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ABSTRACT 

We present new results on the static qq potential from high statistics simulations on 32 4 and 
smaller lattices, using the standard Wilson action at (3 = 6.0, 6.4, and 6.8 on the Connection 
Machine CM-2. Within our statistical errors (~ 1%) we do not observe any finite size effects 
affecting the potential values, on varying the spatial lattice extent from 0.9/m up to 3.3/m. 
We are able to see and quantify the running of the coupling from the Coulomb behaviour of 
the interquark force. From this we extract the ratio y/a/Ai,. We demonstrate that scaling 
violations on the string tension can be considerably reduced by introducing effective coupling 
schemes, which allow for a safe extrapolation of to its continuum value. Both methods yield 
consistent values for A: A.jj^ = 0.558^'^ x y/a = 246_! j MeV. At the highest energy scale 
attainable to us we find a(5GeV) = 0.150(3). 
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1 Introduction 



The experimental determination of the running coupling constant of QCD has reached a 
reasonable degree of accuracy [|l| after two decades of research effort. This has stimulated 
considerable attention to compute this quantity from first principles, by use of lattice 
methods 0, [3|, |J. The lattice approach to the problem of matching perturbative and 
nonperturbative aspects of QCD is notoriously difficult because of the requirement of a 
high energy resolution. Nevertheless, computer experiments in pure SU(2) and SU(3) 
gauge theory have reached a precision that allows to ask rather detailed questions about 
the static quark- ant iquark potential. The size of the available lattices (48 3 x 56, in SU(2) 
gauge theory ||) enables one to decrease the lattice spacing a into a regime where one 
can make contact to predictions of continuum perturbation theory. This has been done 
for the case of SU(2) by a study of the Coulomb behaviour of the interquark force in 
ref. ||. In the case of SU(3), a lattice spacing of a -1 = 3.6 GeV was achieved so far || 
on a 32 4 lattice at (3 = 6.4. This resolution is about the treshold for running coupling 
effects to become visible. 

In this paper we want to present a detailed investigation of the running coupling in 
SU(3) gauge theory, by further reducing the lattice spacing to a -1 = 6.0 GeV. Within 
our analysis of the small distance regime, we will use a parameterization incorporating 
lattice effects. Being limited to lattice sizes up to 32 4 , we have to make sure that our 
results are not spoiled by finite size effects. For this reason we have worked on a variety 
of lattices, at each value of (3. 

Once the running coupling has been extracted, we will be able to compare to pertur- 
bative predictions and estimate a value for the corresponding parameter. We will see 
that this value is consistent with A^, as obtained from the string tension (by the use of 
the two- loop /9-function |7|]), after an extrapolation to a = 0. In order to substantiate this 
result, we will improve on scaling violations (as expressed in the strong /5-dependence 
of A^) by replacing the bare coupling with suitable "effective" couplings |8], [| [II], O , 
measured on the lattice from the average plaquette. In this case, we will find nearly 
asymptotic scaling for (3 > 6.0. The extrapolation to the continuum yields an estimate 
for Al which is consistent within smaller errors with the value obtained from the running 
coupling. 
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2 Methods 



2.1 Sampling 

In order to maintain an appropriate stochastic movement of the gauge system through 
phase space with increasing /3, we have combined one Cabibbo-Marinari pseudo-heatbath- 



sweep |L2] over the three diagonal SU(2) subgroups with four (nine) successive overtax- 



ation sweeps [|13|] for (3 = 6.4(6.8). We reach an acceptance rate of 99.5% for an overre- 
laxation link update. For the heatbath we use the algorithm proposed by Kennedy and 
Pendleton f|TJ| which has a high acceptance rate and can thus be efficiently implemented 
on a SIMD machine. We can afford iterating the algorithm until all link variables are 
changed. On our local 8K CM-2 we need 9.2/zsec for an overrelaxation link update and 
11.5/^sec for a single Cabibbo-Marinari link update. This performance was achieved af- 
ter rewriting the SU (3) matrix multiply routines in assembler language. Measurements 
were started after 2000 - 10000 thermalization sweeps. 



2.2 Smoothing Operators 

In lattice gauge theory physical quantities of interest like masses, potential values, and 
matrix elements are related to asymptotic properties of exponentially decreasing corre- 
lation functions in Euclidean time, and therefore prone to be drowned in noise. So one 
is forced to improve operators in order to reach the desired asymptotic behaviour for the 
small T region. We will shortly describe our particular improvement technique [[J. 

We start from the relation between Wilson loops, W(R, T), and the (ground state) 
potential V(R) 

W(R,T) = C(R)exj){-TV(R)} + excited state contributions. (1) 

Our aim is to enhance — for each value of R — the corresponding ground state overlap 
C(R). Since the ground state wave function is expected to be smooth on an ultraviolet 
scale we concentrate on reducing noise by applying a local smoothing procedure on the 
spatial links: consider a spatial link variable Ui(n), and the sum of the four spatial staples 
IL(n) connected to it: 

H<(n)= J2 Uji^Uiin + jpJin + i). (2) 

3=±1,...,3 

We apply a gauge covariant, iterative smoothing algorithm which replaces (in the same 
even/odd ordering as the Metropolis update) Ui{n) by U-(n) minimizing the local spatial 
action Si(n) = — Re Tr{£/j(ra)nJ(ra)}, which is qualitatively a measure for the roughness 
of the gauge field. This is very similar to lattice cooling techniques already invented 
by previous authors JTR |HJ except that we are cooling only within time-slices and thus 
not affecting the transfer matrix. Alternatively, this algorithm may be interpreted as 
substituting Ui(n) by V (11$ (n)) where V denotes the projection operator onto the nearest 
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SU(3) matrix. In this sense it is a variant of the APE recursive blocking scheme [17] with 
the coefficient of the straight link set to zero, but with even/odd updating. The latter 
feature renders the algorithm less memory consuming and seems to improve convergence. 
Contributions from excited states become increasingly suppressed, as we repeat this 
procedure. After 30(45) such smoothing steps at (3 — 6.4(6.8) we reach values for the 
overlap C(R) of 95(80)% for small (large) spatial separations R. 



2.3 Extraction of Potential Values 

For the extraction of the potential from the Wilson loop data we proceed essentially as 
described in ref. ||, with a slight modification that helps to carry out a straightforward 
error analysis. Instead of fitting the Wilson loops to the dependence 

W S (T- C(R), V(R)) := C(R) exp {-V(R)t} (3) 

for T > T min with some reasonable cutoff T min we take the local mass 

V T (R)=J W, j* r -> ) (4) 
\w(R,T min +l)} 

as an estimator for the potential V(R). By using this explicit formula for the calculation 
of V(R), we are able to propagate the covariance matrix between Wilson loops to a 
covariance matrix for the potential values. This allows one to separate the determination 
of potential parameters from the measurement of the potential itself, helping to decrease 
the degrees of freedom and promoting stability within the fitting procedure. Note that 
the value of V(R) = Vr min (R) does not differ appreciably from the result of a fit to 
eq. (D) because the latter is anyhow dominated by the lowest two T data due to their 
small relative errors. 

The optimization of the overlap C(R) proceeds as described in ref. ||: The parame- 
ters C(R) and V(R) are fitted for different T min to the Wilson loop data separately for 
each smoothing step (and R) according to eq. @ by minimizing 



X %(C(R),V(R))= £ (W^T,) -W^C^MR))) (C™)-^ (5) 



x 



(W(R,T 2 ) -W & (T 2 ;C(R),V(R 



C^? 2 denotes the covariance matrix which is estimated to be 



1 N 

= N(N _ 1) T,(W 1 (R 1 ,T 1 ) - W{R X ,T X )) (WA,T 2 ) - W(R 2 ,T 2 )) . (6) 



i=i 



We have divided the timeseries of Wilson loops into iV successive subsets of given length 
n. Wi(R,T) stands for the average of the respective Wilson loop over the ith subset, n 
should be chosen such that r n N, in order to cope with the autocorrelation time 
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r. Afterwards for each value of R the smoothing step with highest ground state overlap 
C(R) is selected from the fits with reasonable y 2 . 

In a second step stability of local masses Vr{R) (eq. (f|)) against variation of T 
is checked, and T min (R) is determined as the T value (plus one) from which onwards 
stability within errors is observed. For large R values we find T min = 4. For simplicity 
we chose the same value for small R. 

As promised, we are now able to propagate the covariance matrix between different 
Wilson loops C^ 2 to a covariance matrix between the potential values Cy lFt2 , by using 
the quadratic approximation! 

c m 2 = T dV (Ri) c&A dV(R 2 ) 
V & 2 dW(R u Ti) TlTa dW(R 2 , T 2 ) 

T{k\),T[R 2 ) T( J R 1 )+l,T(i? 2 )+l 



W(Rt, T(Ri))W(R 2 , T(R 2 )) W(R l ,T(R 1 ) + 1)W(R 2 , T(R 2 ) + 1) 

riRiR2 /~~tRiR2 
T(/l 1 )+l,T(.R 2 ) T(^ 1 ),T(R 2 )+1 



W{Rt, T{R X ) + l)W{R 2 , T{R 2 )) W(R 1 ,T(R 1 ))W(R 2 , T(R 2 ) + 1) 

where T(R) is used as an abbreviation for T m i n (R). With this covariance matrix we 
are able to fit the potential data to various parameterizations, incorporating all possible 
correlations between different operators measured on individual configurations as well as 
correlation effects within the Monte Carlo time series of configurations. 



2.4 Measurements 

The lattice parameters used for the simulations^ are collected in table |l| which includes 
quotations of 32 4 lattices at (3 — 6.0, and (3 = 6.4, as well as a 24 3 x 32 lattice at (3 — 6.4 
that have been simulated recently ||, and are reanalysed in the present investigation. 
The spatial extent of the lattices at /3 — 6.4 ranges from aLs = 0.87 fm to 1.74 fin. At 
(3 = 6.8 lattice volumes of (0.52 fm) 3 and (1.05 fm) 3 have been realised. The resolution 
a -1 is varied from 1.9 GeV to 6.0 GeV. 

Smoothened on- and off-axis Wilson loops were measured every 100 sweeps (every 
50 sweeps for (3 = 6.0). Up to N max = 30(45) smoothing steps were performed at 
(3 = 6.0,6.4(6.8). The following spatial separations were realized: R = Mei with e« = 
(1, 0, 0), (1, 1, 0), (2, 1, 0), (1, 1, 1), (2, 1, 1), (2, 2, 1). M was increased up to L s /2 for i = 
1,2,4, and up to L^/4 for the remaining directions. Altogether this yields 72 different 

3 In order to check the validity of this approximation, we have moreover carried out a bootstrap 
analysis of our data on the 32 4 lattices. (This method is also shortly described in the appendix 
of ref. ]19[.) The resulting errors (and biased values) are almost identical with the results of our 
approximation, but the bootstrap method alone does not deliver reliable % 2 values (incorporating the 
correlation effects). 

4 Note that we have adapted the physical scales from yfo = 420 MeV to ^fa = 440 MeV. 
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Table 1: The simulated lattices. Physical units correspond to the choice y/a = 440MeV 
for the string tension. Errors ignore the experimental uncertainty within the value of 
the string tension. 





= 6.0 


= 6.4 


= 6.8 


V = L s x L T 


32 4 


16 3 x 32 


24 3 x 32 32 3 x 16 


32 4 


16 3 x 64 32 4 


a/fm 


0.101 (2) 


0.0544 (5) 


0.0327 (5) 


a-'/GeV 


1.94 (5) 


3.62 (4) 


6.02 (10) 


aLs/fm 


3.25 (8) 


0.87 (1) 1.31 (1) 


1.74 (2) 


0.52 (1) 


1.05 (2) 


{aL T y l /MeV 


61 (1) 


113 (1) 


226 (2) 


113 (1) 


94 (2) 


188 (3) 


Total # of sweeps 


6100 


11900 


22000 


10000 


8900 


20400 


15900 


Thermalization phase 


1000 


2000 


2100 


2100 


2500 


10000 


5000 


# of measurements 


102 


100 


200 


80 


65 


105 


110 



separations R on the 32 3 x Lt lattices. The time separations T = 1, 2, . . . , 10 were used. 
Thus the total number of operators measured on one configuration (V = 32 3 x Lt) is 
72 x 10 x N max . 

The potential values at f3 — 6.0, and (3 = 6.4 have been listed in our previous 
publication ||. For convenience of the reader we collect the corresponding values for 
(3 = 6.8 in the appendix. 



3 Results 



3.1 gg-Potential 

We connect our investigation to the recent SU(2) analysis by Chris Michael 0, and 
start from his ansatz: 

V(R) = V + KR-e + I 4ttG l (R)\ + -L. (8) 

The lattice propagator for the one gluon exchange 150] 

- _ r* d 3 k cosjkR) 

has been calculated numerically. The parameter I is expected to be in the range < I < 1 
and controls the violation of rotational symmetry on the lattice (within this ansatz). 
The term f /R 2 mocks deviations from a pure Coulomb behaviour and is expected to be 
positive to the extent that asymptotic freedom becomes visible in the effective Coulomb 
term -(e - f/R)/R. 

A test of the ansatz eq. (|8|) implies that the "corrected" data V(R) = V(R) +SV(R) 
with 

5V(R) = el (A-nG L (R) - 1/r) (10) 
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Table 2: Fit results. Since the parameter values on the largest lattices are most precise, 
we refrain from citing results gained on smaller volumes as long as they are compatible 
with the stated numbers. For the 16 3 x 64 lattice at (3 — 6.8 this is not the case. 
Therefore we have listed both the standard fit result, and the parameter values with the 
string tension constrained to its 32 4 value. 





(5 = 6.0 


P = 6.4 


P = 6.8 


Vol. 


32 4 


32 4 


32 4 


16 3 x 64 


K 


0.0513 (25) 


0.01475 (29) 


0.00533 (18) 


0.00545(27) 


0.00533 


e 


0.275 (28) 


0.315 (15) 


0.311 (10) 


0.269 (22) 


0.274 (18) 


V 


0.636 (10) 


0.6013 (37) 


0.5485 (24) 


0.5412 (37) 


0.5426 (34) 


I 


0.64 (12) 


0.564 (55) 


0.558 (35) 


0.725 (87) 


0.710 (120) 


f 


0.041 (58) 


0.075 (18) 


0.094 (13) 


0.037 (26) 


0.043 (25) 




2 


V3 


^3 


V3 


V3 


x 1 

N nF 


0.816 


0.953 


0.937 


0.989 


0.754 



are independent of the direction of R. The global situation is depicted for the 32 4 lattice 
at (3 — 6.4 in figure 1 where the corrected data points are plotted together with the 
interpolating fit V(R) = Vq + KR — e/R + f /R 2 , with fit parameters Vo,K, e, and / 
as given in table 0. Our potential fits yield x 2 ' /Ndf < 1 as long as the first two! data 
points are excluded. The stability of the string tension result with respect to cuts in R 
is displayed in figure 2 (for (3 = 6.4, and 6.8). 

For (3 > 6.4 the Coulomb coefficients e are definitely different from the value 7r/12 « 



0.262 predicted by the string vibrating picture [21] for large qq separations. The self 
energy contribution V follows the leading order expectation V oc 1/(3. We emphasize 
that for all (3 values the parameter / is established to be positive as expected. In fact, 
this parameter tends to increase with f3, weakening the Coulomb coupling for small 
distances. 

A more sensitive representation of the scatter of the data points around the inter- 
polating fit curve (obtained on the 32 4 lattice) is shown in figure 3 (for (3 = 6.4). Note 
that the deviations are within a 1% band for the largest volume, once the first two data 
points are excluded. Decreasing the lattice spatially or in the time direction by a factor 
of two leaves the data points compatible with the interpolating curve, i.e. the finite size 
effects (FSE) are below our statistical accuracy. Nevertheless it pays to work on a 32 4 
lattice since the larger possible qq separations increase the lever arm needed to fix the 
long distance part of the potential. 

At (3 = 6.8 we find indications of FSE by comparing results from the small lattice 
and the 32 4 lattice. As the string tension appears not to suffer from these effects, we 
have fixed its value to that measured on the larger lattice in order to study FSE on the 
remaining parameters more directly. The largest FSE occurs for the lattice correction 



5 Three for (3 = 6.0. 
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parameter /. This may be due to the low momentum cutoff that starts to become visible 
on the scale of a few lattice spacings. By chosing the form of the one gluon exchange 
(eq. (§)), we have neglected this cutoff in the integral bounds. 

We concentrate our interest here on short distance physics where the linear term is 
not yet dominating the potential. In the case of ft = 6.4 the latter happens at R ~ 5. 
From figure 3 we conclude that reliable results can be extracted for a lattice as small as 
16 3 for this /3-value. In physical units this corresponds to a 27 3 lattice at ft = 6.8. So a 
volume of 32 3 (or even smaller) appears to be sufficiently large for our purpose. 

A synopsis of data for ft = 6.0, 6.4, and 6.8, in physical units, is displayed in fig- 
ure 4 with logarithmic ordinate ranging from 0.03 fm up to 1.9 fm. The three data sets 
collapse to a universal potential. The two curves correspond to a linear plus Coulomb 
parameterization, with the string tension a = Ka~ 2 = (AAOMeV) 2 , and the strength of 
the Coulomb term determined by our fit to the ft = 6.4 data (e = 0.315, full curve), 
and fixed to the Luscher- value (e = 7r/12, dashed curve), respectively. The plot demon- 
strates the incompatibility of the data points with a pure Coulomb behaviour for short 
distances, and the necessity of additional terms like / /R 2 . 



3.2 Running Coupling 

Our lattice analysis for the running coupling a qq (R) closely follows the procedure sug- 
gested in ref. ||. We start from the symmetric discretization in terms of the force F 

a qq (R) = - 3 - Rl R 2 F(R) = ^R lRa YMzmil . (n) 

with R= (Ri + R 2 )/2. We take the corrected potential V(Ri) = V{R\) + 5V(Ri) with 
SV(Ri) as given in eq. (]T0|) . Unlike ref. ||, however, we use all possible combinations 
R U R 2 with \R X - R 2 \ < 1.5. 

The resulting data points are contained in figure 5a. In order to exhibit both the 
global behaviour, and the perturbative region (R — > 0) we decided to use a logarithmic 
ordinate (in units of a -1 / 2 ). The latter region is expanded in the inset. We omitted all 
values with errors Aa qq (R) > a qq -(R)/3 in order not to clutter the graph. In addition to 
the statistical error of the force F(R) we allow for a systematic error 

^F syst (R) = ((^] 2 + (—) 2 ) \5F(R)\ (12) 



with 6F(R) = ^^Zr^ ■ &F syst is typically of the order of 10% of the lattice correc- 
tion 5F(R). 

Now we can proceed to analyse our a^-data in terms of the continuum large momen- 
tum expectation for the running coupling: 

oc qq {R) = (b \n(RaA R y 2 + b 1 /b ln\n(RaA R )- 2 )~ 1 , (13) 
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with 

11 No M(N C \ 2 . . 

being the first two coefficients of the weak coupling expansion of the SU(Nc) Callan- 
Symanzik /3-function (eq. (ESj) below). In order to extract Ar we base our fits exclusively 
on data points at @ — 6.8 with Ri, R 2 > \/S on the r.h.s. of eq. (pTTf) . This is done in 
order to avoid the danger of "pollution" from discretization errors. 

We now ask the question, within which R region our data are compatible — if at all! 



- with the asymptotic behaviour of eq. (|T3|) . We find that as long as RyK < 0.173 our 
fits yield results with reasonable y 2 /Np>p. This upper limit in R corresponds to 2.5 GeV. 
Fitting the j3 = 6.8 data over this region we obtain 

A R = (0.562 ± 0.020 ± 0.010)^ « (247 ± 10)MeV. (15) 

The first error stems from the fit just described, while the second relates to the statistical 
uncertainty of the string tension within our lattice analysis. The corresponding fit curve 
with error bands is plotted in figure 5. As the data appear to osculate the asymptotic 
curve one finds a systematic dependence on the R cut: A# tends to be larger if more 
(low energy) data points are included and vice versa. In this sense one might consider 
our value as an upper limit. 



Exploiting the relation Ar = 30.19Ai we get: 

A£ = (18.6 ± 0.7 ± 0.3) x 10~ 3 ^ « (8.19 ± 0.33)MeV. (16) 
This corresponds to the ratio 

= 53.7 ±2.1. (17) 



a 
At 



In figure 5b we have plotted a versus the energy. At the largest realized energy scale we 
find a qq (5GeV) « 0.150(3). 

Returning to the global structure of the data displayed in figure 5 we make three 
observations: 1. The small R contributions (circles, and triangles) follow very neatly the 
asymptotic perturbative prediction eq. (|13D , indicating very little discretization effects. 
2. Over the whole R range the data sets for (3 = 6.4, and (3 = 6.8 coincide very nicely, 
giving evidence for scaling. 3. The deviations of the data from the asymptotic behaviour 
remain fairly small up to q « lGeV or a q q ~ 0.4. 

We conclude that lattice simulations can indeed make contact to the perturbative 
regime. Moreover, it is very satisfying to observe that the 2-loop-formula describes the 
lattice data down to a scale as small as 1 GeV — at least in the quenched approximation 
of QCD. One would expect that the situation in full QCD is fairly similar, concerning this 
property. In the infrared regime (q < y/a) the differences between both theories will be 
considerable: Because of the linear confining potential our expectation for the pure gauge 
sector is ct q q(q) oc l/q 2 . This has to be confronted with the expression oi qq -(q) oc e~~^ g /q 2 
for QCD with fermionic degrees of freedom where \x stands for the screening mass. 
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Table 3: The lattice spacing a, and cutoff parameters Al calculated from the 
2-loop-expansion eq. in units of the string tension a. Al is obtained by insert- 
ing the bare lattice coupling. For ' the j3 E ' effective couplings were used. A naive 
linear extrapolation to a = leads to the results displayed in the second last row. Log- 
arithmic extrapolations yield the values in the last row. 



p 


ay/a 


v^/Al 






5.7 


0.4099 (24) 


124.7 (0.7) 


63.3 (0.4) 


55.7 (0.3) 


5.8 


0.3302 (30) 


112.4 (1.0) 


63.0 (0.6) 


55.6 (0.5) 


5.9 


0.2702 (37) 


102.9 (1.4) 


61.2 (0.8) 


54.3 (0.7) 


6.0 


0.2265 (55) 


96.5 (2.3) 


60.0 (1.5) 


53.4 (1.3) 


6.2 


0.1619 (19) 


86.4 (1.0) 


56.9 (0.7) 


50.8 (0.6) 


6.4 


0.1215 (12) 


81.3 (0.8) 


55.7 (0.5) 


50.0 (0.5) 


6.8 


0.0730 (12) 


76.9 (1.3) 


55.7 (0.9) 


50.4 (0.8) 


oo 


lin. 





63.6 (2.4) 


53.1 (1.6) 


48.3 (1.4) 


log. 





54±ig 


53.2 J g 


49. lt^ 



3.3 Scaling 

Normally one speaks of asymptotic scaling when the ratio \fa / Al remains constant 
on varying f3 where 

^ =1 a^{-2^)^' k (18 > 

(with g 2 = 2N c //3) denotes the integrated two- loop /3-function (eq. ( P5| ) below). In 
table [3] we have compiled our new results on the string tension together with previous 
results from refs. || ^3|. As can be Seen, we are still far away from the asymptotic 
scaling region up to (3 = 6.8. 

We attempt to extrapolate A^ 1 to the continuum limit by the use of a parame- 
terization that takes into account the leading order expectation for scaling violations 
0(1/ In a): 

A7» = A- L \0) + g- (19) 

yja m(Day/a) 

We find the data compatible with this logarithmic behaviour, with D pa 1-2, and C ~ 
20-80. The fit parameters are not particularly stable with respect to a variation of 
the number of data points. The bandwith of extrapolations to the continuum limit is 
illustrated in figure 6a where we have plotted the extreme cases of a fit to our four 
low a data points, and all seven data points (open circles). If we average the values 
obtained from these fits, and take the upmost and the lowest possible numbers as error 
bandwidth, we estimate the asymptotic value to be a/cxA^O) -1 = 54jl 15 (full circle). We 
would like to mention that a naive linear extrapolation to the continuum limit yields the 
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Table 4: The average plaquette action (Su), measured on large lattice volumes. The 
values for (3 < 5.9 are taken from the collection in ref. [[□]] while the other numbers are 
our new results, obtained on 32 4 lattices, and one 24 3 x 32 lattice (f3 = 6.2). 



p 


(Su) 


5.7 


0.45100 (80) 


5.8 


0.43236 (5) 


5.9 


0.41825 (6) 


6.0 


0.406262(17) 


6.2 


0.386353 (8) 


6.4 


0.369353 (5) 


6.8 


0.340782 (5) 



value v /aAi(0) 1 = 63.6(2.4) with (obviously) underestimated error. We take this as a 
warning for purely phenomenological continuum extrapolations. 

In view of the uncertainty of the above number it would be highly desirable to improve 
the situation by developing a scheme within which the a dependence of Ax, (a) is reduced. 
Parisi suggested many years ago a more "natural" expansion parameter g% ||, based on 



a mean field argument. His scheme was elaborated in refs. P, [TOl 11 1. It works as follows: 



Let c n be the coefficients of the weak coupling expansion of the average plaquette 



oo 



E Cn9 2n - (20) 



n=l 



The idea, now, is to introduce an effective coupling in terms of the Monte Carlo generated 
average plaquette 



9% = ^ (21) 

Cl 

Ci Ci 

for which the first order expansion is exact. The hope is that the nonperturbative (or 
higher order perturbative) contributions that are resummed in the effective coupling gg 
may compensate high order terms in the /5-function which are responsible for the scaling 
violations. Support for this expectation comes from the observed scaling of ratios of 
physical quantities (figures 4, 5) within the same f3 region. 

The coefficients c\, and C2 have been calculated previously |24], and an unpublished 



value for C3, obtained by H. Panagopoulos, has been cited in ref. [II]. The numerical 
values are: 

ci = (Nq — 1)/ (8Nc) (22) 
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c 2 = (N 2 - 1)(0.0204277- l/(32A^))/4 (23) 
c 3 = (JV^-l)JV c (0.0066599-0.020411/JV^ + 0.0343399/iV^)/6. (24) 

The plaquette values needed for the conversion into the effective coupling schemes are 
collected in table |j. The numbers for j3 < 5.9 were taken from the collection in ref. [11 
Starting from the expansion 



M = = - E Kg*** (25) 



of the /3-function, one rewrites 



dg E dg g dg 



dlna dlnagE dg 
-b o9 % - hg E - b 2 g 7 E + Ubo (2 (^) 2 - f) - 2hf) g\ + 0(g%) 



The first two terms in this weak coupling expansion remain unchanged under the 
substitution. Therefore, an integration again leads to eq. (|T8|), but with a redefined 
integration constant 

A E = A L exp(-^-) ^ 2.0756A L (for SU(3)). (27) 
\2cib / 

This factor is due to a shift of the effective (3 by a constant in the continuum limit: 
g% 2 = g~ 2 — c 2 /ci + 0(g 2 ). In the following we will refer to this scheme as the (3^ 
scheme. As one can see from figure 6a (open squares), and table |3| this kind of (numerical) 
resummation of the asymptotic series eq. (p0| ) leads to considerably reduced logarithmic 
corrections (C ~ 2.5). 

As an additional check of this improvement technique we consider in the following 

(2) 

an "alternative" effective coupling scheme j3 E . Our idea is to introduce a coupling g 2 
by inverting the relation 

(S D ) = c ig 2 2 + c 2 g A 2 . (28) 

This amounts to truncating the weak coupling expansion eq. ( p0|) after the second term!. 
A short calculation yields: 

13(92) = -bo9 3 2 - hgl - b 2 g\ - 3b % 7 2 + 0(g 9 2 ). (29) 

Cl 

Because of g 2 2 = g~ 2 + 0(g 2 ) the integration constant Ax, remains unchanged in respect 
to the original bare coupling scheme. 

6 One can generalise this scheme by truncating in higher orders n. This is of little interest, however 
(unless one is interested in numerical studies of the impact of a particular higher loop contribution on 
the observed scaling violations), since the /^-function has only been calculated up to 0(g 5 ). Moreover, 
one would retrieve the bare coupling scheme at n sufficiently large. 
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If we compare the third order terms of the two effective schemes (eqs. ([26 ,29)), we 
find explicitely: 

M = (3(9e(9)) + 5.3 x 10-y + 0(g 9 ) = (3(g 2 (g)) + 4.02 x lO'V + 0(g 9 ). (30) 

This means that the correction of the /3-function through the 3-loop-contribution is much 
larger for the (3^ than for the (3^ schemel. Nevertheless, at least within the investigated 
(3 region, the qualitative behaviour of both schemes is the same as can be seen in figure 6a. 
For the (3^ scheme the correction coefficient (C ~ 0.9) of the continuum extrapolation 
eq. ([19]) is even smaller than for the (3^ scheme. In figure 6a we have included the 
estimates for the asymptotic A^ 1 values (and the A^ 1 from the running coupling) as full 
symbols. 

The extrapolated values for both effective schemes are, respectively: 

= 53.2^A« (31) 
= 49.lj£jjA?. (32) 

Averaging these numbers that carry asymmetric errors leads to \fa = 50.8^4;gAf . This 
result is in nice agreement with the ratio extracted from the running coupling (^fa = 
53.7(2.1)A2, eq. (|i~7|)). Using this additional information, we obtain: 

V^=51.6l?;£A L . (33) 

This result may be converted into any continuum renormalization scheme like the min- 



imal subtraction (MS) scheme. By exploiting the relation A^r = 28.81A^ [p5[| , we 
get: 

^ = 0.558lffi. (34) 
v ° 

Let us finally comment that the two approaches presented in this paper for the 
determination of the QCD scale parameter A, namely to analyze (a) g 2 (AaR), and (b) 
its inverse Aa(g 2 ) in terms of the two-loop predictions eqs. (p!3|JT8|) , are complementary 
and supportive to each other because higher order corrections to methods (a) and (b) 
are anticorrelated. In our running coupling (string tension) analysis we observe the 
"effective" A^ to decrease (increase) with the energy scale. Since the central value of 
our "upper limit" A£ is smaller than that of our "lower limit" Af we are in the position 
to state relatively small errors for Aj^. 

In figure 6b we have plotted the A^ 1 data versus (3 in order to visualize the slow 
approach of the bare coupling data towards the asymptotic value, and the improvement 
achieved by the use of effective couplings. 



7 Note that the difference between the /3-functions for both effective schemes is independent of C3 to 
this order. 
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4 Discussion 



We have demonstrated that medium size computer experiments are able to determine the 
A-parameter of SU (3) Yang-Mills theory within a reasonable accuracy (that can compete 
with QCD experiments). For this result, it has been important to study both infrared, 
and ultraviolett aspects in order to verify the reliability of the continuum extrapolation. 
We might say that we have been lucky to get hold of asymptotia within our means. This 
is due to the discovery that the running coupling constant is well described within this 
theory by the two-loop formula down to a scale of about 1 GeV. 

If nature continues to be nice to us, and the inclusion of dynamical quarks results 
only in a /3-shift of quenched predictions it is possible to predict experimental numbers 
like as(M z ), as explained in ref. 0. Obviously, it is preferable to repeat this study in 
full QCD on the level of TERAFLOPS power. In the meantime, further improvements 
of lattice techniques are of great interest. A promising route has been proposed by 
M. Liischer et. al. [p6 |, and tested on SU(2) Yang-Mills theory. These authors start 



from a volume dependent coupling g(L) which allows them to reach large energies on 
small lattices. 

After completion of this work we received a preprint by S.P. Booth, C. Michael, and 
collaborators |27| that contains a running coupling study for 577(3) gauge theory up to 
(3 = 6.5. Their results are fully consistent with ours. 



Acknowledgements. We are grateful to Deutsche Forschungsgemeinschaft for the sup- 
port given to our CM-2 project. We thank Peer Ueberholz and Randy Flesch for their 
kind support. One of the authors (G.B.) would like to thank Chris Michael, Edwin Laer- 
mann, Rainer Sommer, and Jochen Fingberg for helpful discussions about data analysis, 
and the different effective coupling schemes. 
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Appendix 



A Potential Values 



In this appendix we are stating the potential values measured on a 32 4 lattice at f3 — 6.8. 
The corresponding numbers for the other j3- values can be found in ref . || . The on- and 
off-axis paths are numbered in the following way: 



Path # 


1 


2 


3 


4 


5 


6 


Path (X, Y, Z) 


(1,0,0) 


(1,1,0) 


(2,1,0) 


(1,1,1) 


(2,1,1) 


(2,2,1) 


Elementary distance M 


1 


1.41 


2.24 


1.73 


2.45 


3 



The results for the potential V(R) (in lattice units), as well as for the "corrected" 
V(R), and the corresponding ground state overlaps C(R) are collected in table 5. The 
data is plotted (among the other curves) in figure 4. 

Table 5: The potential values V(R) (in lattice units a -1 ), "corrected" values V(R), and 
groundstate overlaps C(R) for (3 = 6.8, V = 32 4 . 



R 


Path 


V(R) 


V(R) 


C(R) 


1.00 


1 


0.3107 (6) 


0.3210 (10) 


0.950 (3) 


1.41 


2 


0.3855 (11) 


0.3794 (12) 


0.951 (4) 


1.73 


4 


0.4188 (19) 


0.4098 (20) 


0.946 (8) 


2.00 


1 


0.4236 (14) 


0.4266 (14) 


0.929 (5) 


2.24 


3 


0.4428 (13) 


0.4397 (14) 


0.934 (5) 


2.45 


5 


0.4559 (15) 


0.4509 (15) 


0.936 (6) 


2.83 


2 


0.4696 (20) 


0.4656 (20) 


0.923 (8) 


3.00 


1 


0.4725 (14) 


0.4709 (14) 


0.931 (6) 


3.00 


6 


0.4751 (18) 


0.4705 (19) 


0.924 (7) 


3.46 


4 


0.4906 (31) 


0.4861 (31) 


0.923 (12) 


4.00 


1 


0.5000 (18) 


0.4970 (19) 


0.916 (7) 


4.24 


2 


0.5079 (23) 


0.5039 (23) 


0.939 (9) 


4.47 


3 


0.5105 (22) 


0.5068 (22) 


0.916 (9) 
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Table 5, continued 



± V 


Path 

X d 111 


V(R) 

\ / 


V(R) 

\ / 


CAE) 


A 90 


d 


51 78 (28^ 

VJ . d X 1 (J IZj(JJ 


9 51 39 (98"l 

VJ.5jXO<J 1 Z(J 1 


n 91 3 fi 1 

U. Jld 1 11 1 


5 nn 

d .VJVJ 


1 


n 51Q3 


n 51 59 (90) 

VJ.dXdC/ I ZU 1 


n 994 (R\ 

VJ. C/ZjT: 1 O 1 


5 2D 

• J . Zj VJ 


4 

T: 


n 523n (32^ 


n 5199 f32"l 

VJ.tJXiyVJ 1 iJZ 1 


9 999 3^1 

VJ . U £iU 1 Xd ( 


5 66 

d . VJVJ 


9 

Zj 


n 531 2 (29"l 

U.UUlZ I Zj <J J 


n 5273 f3n"i 

VJ.tJZjJO IdVJI 


n 99n o 9) 


6 nn 


1 

X 


n 5325 f25"l 


n 5289 (95"l 

Kj.xJ^ikjU IZjdf 


n 997 0) 

VJ . <J VJ f I XVJ f 


6 nn 


6 

VJ 


n 5357 f29"l 

VJ.tJOtJ 1 1 Zj<J J 


n 531 7 (30^ 

VJ.ddX 1 IdVJf 


n 91 8 9) 

VJ.«JX(J IXZjJ 


6 71 


d 


n 5421 (97) 

VJ.dlZjX 1 Zj 1 J 


n 5383 (97) 

VJ.tJO(JO 1 Z 1 J 


n 9i 7 ('\'\) 

VJ . <J X 1 I X X ( 


6 93 


4 

T: 


n 5469 (42"l 

VJ.dTiVJ<J I TiZj i 


n 543n r42"i 

VJ.dTidVJ 1 1Z 1 


9 91 6 fl 6^1 

VJ.tVXVJ IXVJI 


7 nn 

1 .VJVJ 


1 

X 


n 5463 (97) 

VJ .dTiUd 1^1 J 


n 5496 (97) 

U .tJt:ZjVJ 1 Zj 1 J 


n 991 (\ 1 

VJ . i7Zj X 111 1 


7 n7 

I .VJ I 


o 


n 5474 (If,) 

VJ .dTi 1 t: 1 OU 1 


n 5436 (17) 

VJ .dTidVj Id 1 1 


n 998 (\ 5^1 

VJ. C/ZjO 1 Xd ) 


7 35 


d 


n 55n4 f39"l 


n 5466 f39"l 

VJ.dTiUU I dz 1 


n 993 3"l 

VJ. C/Zjd 1 Xd 1 


8 nn 

(J . VJVJ 


1 

X 


n 5568 (34"l 

VJ.ddVJw IutI 


n 5531 (34"l 

VJ.dtJOX 1 OtI 


n 91 n fi 3"i 

Vj.iyXVj 1 iu 1 


8 49 


9 


n 5623 (44"l 

\J.xJ\J^ixJ \ tit: 1 


n 5584 (44"l 

VJ.dddTi 1 tit: 1 


n 91 1 fi 7) 

VJ . <J X X 111 J 


8 66 

O . VJVJ 


4 

T: 


n 5644 f47l 

VJ.dVjTiTi 111 J 


n 56n5 (A7) 


n 93n (i q) 

VJ . .JdVj I 1 J I 


8 Q4 

O . C/t: 


d 


n 5663 f37"l 

VJ .dUUd Id 1 1 


n 5695 ('37 N I 

VJ .dUZjd Id 1 1 


n 91 1 n 5"i 

U. C/ X X I Xd ) 


q on 


1 

X 


D i(i71 f^fi^ 

VJ . dVJ 1 X I dVJ J 


\J.O\JOO \0\J J 


Q90 (AA\ 

VJ . i^ZjVJ I 11 J 


q nn 

u .\J\J 


6 

VJ 


n 5651 (37^ 

VJ .dUd x \o i J 


n 561 9 (37"l 

VJ .dU X Zj Id 1 1 


n 91 1 n 5"i 

U. C/ X X I Xd ) 


q sn 

<J . OVJ 


d 


n 5733 (41 

VJ.d ( y^" 1 - 1 - / 


n 5695 (41 

VJ.dVJtyd I t1 J 


9 998 <n) 

VJ . <J VJ(J 1 X VJ 1 


q nn 

<J . <JVJ 


2 

Zj 


n 5745 f48 

VJ.d ( JT "J I1U J 


9 57n7 (48"l 

VJ . d 1 VJ 1 1 t:(j J 


n 995 f1 Q^i 

VJ . JZU 1 U 1 


1 n nn 

-L VJ . VJVJ 


1 

X 


n 5743 f44"l 

VJ.d 1 T-J I tit: J 


n 57n5 (44"i 

VJ.^JiVJ^J \ tit: 1 


n 994 n 7^1 

VJ . <J VJt: 111 ( 


i n 39 


4 

T: 


n 5777 (53^1 

U.d 1 1 1 1 dd f 


n 5739 l'53'l 

yj.^Ji^JU Iddf 


n 995 T91 

VJ . <J VJd 1 Zl 1 


1 1 nn 

1 1 . uu 


1 

X 


n 583n f4Q"l 

\J.OOO\J 1 TiC/ 1 


n 5799 r49"l 

VJ .d 1 C/Zj IIj J 


n 91 3 (1 9^1 

VJ.C/ld 1 Ij 1 


X X . XO 


o 


n 581 8 f4Q"l 

U. JOIO 1 TiC/ 1 


9 578n (AQ) 

VJ .d 1 OU llj 1 


n 993 (90) 

VJ. C/VJd 1 ZjVJ J 


1 1 31 

X X . d X 


9 

Zj 


n 5841 f5n"i 

U.UOtiI 1 dU 1 


n 58n3 r5n"i 

VJ.dOVJd 1 dVJ I 


n 898 f9n"i 

VJ. OC/O 1 ZjVJ J 


1 9 nn 


1 

X 


n 5887 (55^ 

VJ .dOO I \oo J 


n 5849 (55"l 

VJ.dOTiC/ 1 dd 1 


n 895 (91 

\J .OCfO 1 Zj 1 1 


1 9 nn 

XZj. uu 


6 

VJ 


n 5nnn i'55'i 


n 5869 (55"l 


n 9m (99) 

VJ . <J VJ 1 I Z(Z( J 


12.12 


4 


n 5941 f6n) 

u.Ul/t:i I \J*J / 


n 5992 (50) 


n 928 f24) 


12.25 


5 


n.5918 (53) 


n.5879 (53) 


0.912 (21) 


12.73 


2 


n.5962 (64) 


9.5923 (64) 


0.918 (26) 


13.00 


1 


n.5987 (56) 


9.5949 (56) 


0.912 (22) 
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Table 5, continued 



± V 


Path 

X d uJX 


!/(/?! 

\ / 


V(R) 

\ / 




1 3 42 

A.O . TlZj 




n 5qqs (fii 1 

VJ .OU CtO I VJ X 1 


o 5n6n (61 l 

W.*J<JVJW I VJ X J 


o 8n4 r24"i 

VJ.OiJt: I ZjT: 1 


1 3 8fi 


4 

T 


n fi03i f75l 


n 5qq3 (751 

VJ .OctctO \ I O ) 


n 8n5 f2Q"i 

VJ. OcrO I ZjC/ 1 


14 nn 


1 


n fin55 (621 

W.VJWtJtJ IVJZvJ 


n 6ni 7 (fi2i 

VJ.VJVJX ( IVJZj J 


n sqq r24"i 

\J.\JUU I Zjt: 1 


14 14 

It: . 1 T 


9 

Zj 


n fin52 (731 

VJ.VJWtJZj II U 1 


n 6ni 4 (731 

VJ.VJVJXt: II Ul 


n 8n3 r2q"i 

VJ.O.JO IZviyJ 


14 7D 

It:. 1 W 


tj 


n finqfi (701 

W .VJW <JVJ I 1 VJ J 


n 6n58 (70) 

VJ .VJVJUO \ i VJ J 


n 8n5 ('271 

VJ . OiJ ^J 1^1 J 


1 5 nn 


1 

X 


n finn7 (fi8l 


n 6n5q (fi8l 

VJ.VJVJtJiJ I UU 1 


n 8n5 (27) 

vj . OiJ ?j l^i y 


1 5 nn 


fi 


n fii n2 (fiql 

U.U1UZ \VJUt 


n 6n64 (691 

VJ.VJVJVJt: I VJ <J J 


n 8n5 (27) 

VJ . OiJ ?J 1^1 J 


1 5 56 

XtJ . 0\J 


9 

Zj 


n fii 3q (81 ) 


n 61 m (81 l 

VJ.VJXVJX 1 Ul 1 


n qn4 (12) 

VJ.tyVjT: IOZjJ 


1 5 5Q 

J.O .Oct 


4 

T 


n fii fi3 (8i l 

VJ. U X VJO I O X 1 


n fii 25 (81 1 

VJ .U X Z-tJ \ O X 1 


n qi n f32l 

VJ. C/ XVJ I OZj f 


1 5 fi5 


Q 


n fi1 44 (731 

VJ . VJ X Hit: I 1 O 1 


n fii nfi (731 

VJ .VJ XVJVJ \\Oj 


n 8n5 r2qi 


i fi nn 


1 

X 


n fi1 51 (741 

VJ.UXtJX 114 1 


n fii 1 3 f74 1 

VJ.VJXXO NtI 


n 878 (9Q) 


1 fi Q7 


2 

Zj 


n fi24fi (881 

VJ . VJZjtiVJ 1 (JO 1 


n fi2nq (rr) 

VJ.VJ^jVJU IOOI 


n 886 ('341 

VJ.OOVJ 1 Oil 


17 15 


*J 


n fi248 (7R) 

VJ.VJZvTU N U 1 


n fi2i n (7R) 

VJ .VJ 4j ±VJ 1 1 U 1 


n 8n5 f3i 1 

VJ.OiJtJ 1 ul 1 


1 7 32 


4 

T 


n fi258 (Q4-) 

W.VJZvtJO I Jt 1 


n fi22n (q4i 

W.VJZvZvW 1 Jtl 


n 883 ('361 

VJ.OOO lOVJJ 


1 7 8Q 




n fi2qfi fqol 

U.UZjC/U \CtVjj 


n fi258 ('QO^i 

U.UZjUO \CtVj J 


n 88n f35l 

VJ. OOVJ I OO J 


1 8 on 


u 


VJ.VJOXZv 1 OO J 


fi974 ^RR^i 

VJ.VJZvIt: 1 OO 1 


n 885 (341 

\J.(J(JO 1 Ot 1 


1 8 3Q 


9 

Zj 


n fi337 ('QQ^i 

VJ. VJOO 1 I C/C/ 1 


n fi2qq (qql 

VJ.VJZjC/C/ \CtctJ 


n sqq f3ql 

VJ. Octet \Octj 


1 q n5 

-L iJ . VJO 


4 

T 


n fi3q4(i oql 

VJ . Ul) Jt I X VJ <J 1 


n fi357(1 DQ^l 

VJ.VJOtJl I 1UJ 1 


n qnn T43l 

vj.iyvjvj iiui 


1 q fin 


tj 


n fi4n2 (q5l 


n fi3fi4 (o^) 

VJ.VJOVJt: ItJ^JJ 


n 874 (361 

VJ . O I t: \O\Jl 


1 q 80 


2 

Zj 


n fi44on 051 

VJ . VJttVJ I X VJO f 


n fi4n2fi n5"i 

VJ . VJtVJZj I X VJtJ ( 


n 882 f41 1 

VJ . OOZj 111 1 


2n 7q 

ZjVJ . ( <J 


4 

T 


n fi48fi(4 1 7) 

VJ . VJtOVJ I X X I J 


n fi448f1 1 7) 

VJ . VJttO 1111 1 


n 875 f44l 

VJ . O 1 tJ I TT I 


21 nn 

Zj J. .uu 


fi 


n fi4qfiH 081 

VJ. VJtC/VJ I XVJO 1 


n fi458Tl ns^i 

VJ.VJttJOI XVJO 1 


n 87q T421 

U.O I c? v^^' / 


21 21 

Z, X . Z, X 


9 

Zj 


n fi52fif1 1 81 

VJ. UUZjU I X XO 1 


n fi48nn 1 8"i 

VJ.VJtOC/I X XO 1 


n 8m T4fii 

U.OC/ 1 ItU 1 


22 52 

. U Z, 


4 

T 


n fifii on 31 i 

U.UUXUIXOX 1 


n fi573f1 31 

VJ.UtJl Ol lOl 1 


n 88q f5i 1 

\J.O(Dct lull 


22 fi3 

ZvZv .VJO 


9 

Zj 


n fi545(4 281 

VJ .VJO'-kO \ X ZjO 1 


n fi5n8H 28"l 

VJ.VJtJVJOlXZjO 1 


n 847 f48l 

U.Ot: 1 ItU 1 


24 nn 

Zjt: .WW 


fi 


n fifi88(423l 


n fifi5nfi23^ 

VJ . \J \JtJ VJ I 1ZU ( 


n 8fi3 f47l 

VJ.UVJO I t: I 1 


24.25 


4 


n.66q4(142) 


n.6657(142) 


n.862 (53) 


25.q8 


4 


0.6791(151) 


0.6753(151) 


n.866 (57) 


27.71 


4 


n.6qn8(i62) 


0.6871(162) 


n.848 (60) 
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